Low prevalence, quasi-stationarity and power-law distribution in a model of spreading 
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Understanding how contagions (information, infections, etc) are spread on complex networks is 
important both from practical as well as theoretical point of view. Considerable work has been 
done in this regard in the past decade or so. However, most models are limited in their scope and 
as a result only capture general features of spreading phenomena. Here, we propose and study 
a model of spreading which takes into account the strength or quality of contagions as well as 
the local (probabilistic) dynamics occurring at various nodes. Transmission occurs only after the 
quality-based fitness of the contagion has been evaluated by the local agent. The model exhibits 
quality-dependent exponential time scales at early times leading to a slowly evolving quasi-stationary 
state. Low prevalence is seen for a wide range of contagion quality for arbitrary large networks. 
We also investigate the activity of nodes and find a power-law distribution with a robust exponent 
independent of network topology. Our results are consistent with recent empirical observations. 



PACS numbers: 89.75.Hc, 05.70.Ln, 87.23. Ge, 89. 75. Da 



INTRODUCTION 

Spreading, defined broadly, as transmission of conta- 
gions (e.g. information or virus) from one agent to the 
next is a fundamental process with application in a wide 
range of disciplines including physics, epidemiology and 
social sciences. In many cases of spreading phenomena, 
transmission only occurs after some local conditions are 
met. From a theoretical point of view, this requires 
construction of models of contagion dynamics (spread- 
ing) which take into account the content of the conta- 
gion as well as the role of the individual agents receiv- 
ing/evaluating/transmitting contagion. It is therefore 
surprising that such models have not been proposed and 
studied, despite the fact that much attention has been 
given to spreading phenomena, typically in a mean-field 
approximation, in the past decade (lHiJ. In this Arti- 
cle, we propose and study a model of contagion spread- 
ing which takes into account the "quality" of contagion 
as well as the role of the individual local agents in the 
spreading phenomenon. We find many interesting and re- 
alistic features in our model, including quality-dependent 
fast initial spreading followed by extremely slow dynam- 
ics, as well as robust power-law behavior in activities of 
individual agents. Low prevalence is also observed as a 
generic but limiting behavior of our model, e.g. as system 
size diverges. 

It is hardly possible to overemphasize the importance 
of contagion dynamics. Modern telecommunications in- 
cluding emails and text-messages, and more recently so- 
cial networks like Facebook and Twitter have revolu- 
tionized the process of information spreading. Network 
and/or viral marketing [f|, opinion formation and ru- 
mor/innovation spreading [|| 0], recruiting and talent 
searches [8| .disaster response and relief efforts, mobiliz- 
ing masses [§] , as well as epidemic spreading [l|, 0, EH are 
a few examples of how contagion spreading through com- 



plex networks is of vital importance in our modern way 
of life. Scarcity of reliable empirical results, as well as the 
ease to capture the general features of spreading, have led 
many authors to study spreading along the lines of epi- 
demiology P. [lll - flij . In such approaches, one typically 
uses SI, SIR, or SIS models on a complex network where 
the letters in the acronyms refer to the state (Suscepti- 
ble, Infected, or Recovered) of the agents on the network. 
Assuming a mean-field transmission rate A, one typically 
finds a low threshold for spreading where a large part of 
the network is "infected" in a relatively short time, with 
some dependence on network topology 0, [l3l - 15 1 . How- 
ever, such results are in contrast with most real- world 
observations which indicate a low prevalence in a variety 
of spreading phenomena [16| . Here, we introduce the con- 
cept of contagion "quality" in such models and take into 
account the role of individual local agents in evaluating 
transmission condition based on the observed "fitness" of 
such contagion. We employ local probabilistic conditions 
for transmission on two types of most commonly studied 
complex networks. We find that our model can provide 
a more realistic picture of spreading phenomena, leading 
to results consistent with general empirical observations 
of such phenomena. 



A MODEL WITH CONTAGION QUALITY AND 
LOCAL DYNAMICS 



In most commonly studied models of spreading, the 
role of the active (local) agents is often reduced to a 
uniform (mean-field) behavior where upon receiving, the 
local agent transmits the contagion with probability A 
independent of the particular agent. Some important 
(but less- noticed) exceptions may be found here [17H19| . 
Another shortcoming of such general models is the fact 
that the "quality" or "strength" of the contagion is often 
neglected. However, not all contagions are equal. For 
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example, a new virus or a shocking news have a much 
more probability of transmission than a common virus or 
a normal everyday news. In fact, both these shortcom- 
ings of the most frequently studied models are related in 
a sense that in most spreading phenomena, local agents 
"judge" the quality of the transmission received before 
deciding to pass it on. Typically, the contagion interacts 
with the local agent, and if some "fitness" conditions are 
met, transmission occurs. To the best of our knowledge, 
models where local agents use the quality of the conta- 
gion they receive in order to decide whether or not to 
pass it on have not been studied before. In the present 
work, we introduce a quality factor for the contagion be- 
ing spread and thus introduce local dynamics whereby 
various agents act differently based on the quality of the 
contagion they receive. We find that the spreading pro- 
cess is characterized by an initial phase where conta- 
gion spreads exponentially fast with time scales which 
strongly depend on the quality of the contagion being 
spread. This leads to a quasi-stationary state where ex- 
tremely slow evolution eventually leads to a final equilib- 
rium state. We note that such slow spreading occurring 
after the initial exponential growth has been observed in 
various contagion spreading phenomena before 2ol. 21 1. 
However, the cause of such slow spreading has been as- 
sociated with various types of correlations that may ex- 
ist within standard models (2lj . In our model, however, 
both exponential and slow spreading exist with quality- 
dependent time-scales, and this general behavior is shown 
to hold on various network structures regardless of topol- 
ogy. Another important and often neglected feature of 
spreading is the activity of the agents. We therefore in- 
vestigate such quantity and show that for a wide range 
of contagion quality, agent activities exhibits a power-law 
distribution, which again occurs independent of network 
structure. Such power-law behavior has recently been 
seen in real social networks 2(3-23}. 

Accordingly, we introduce a local quantity (called qual- 
ity) Xi for each node i. Here, we assume that the quality 
of each individual is directly proportional to the number 
of its neighbors fej, x\ = ki/k max where k max = max(ki), 
making < Xi < 1 for any given network. We also in- 
troduce a parameter a which characterizes the quality 
of the contagion being spread, where < a < 1. Here, 
we consider a local probabilistic rule for acceptance (i.e. 
transmission) of contagion based on the perceived fitness 
of the incoming contagion, a, which is defined as follows: 
for each node i which receives the contagion a, with prob- 
ability f(ki,a) — exp(—(xi — a) 2 /2a 2 ) accepts the con- 
tagion and with probability 1 — f(ki,a) rejects it. This 
is essentially our fitness criterion. Note that if the con- 
tagion is accepted, the individual keeps that contagion 
forever and passes it on to all its ki neighbors, otherwise 
the individual rejects the contagion. We note that other 
forms of stochastic as well as deterministic a-dependent 
local rules may be considered [24j |. Here, we consider such 



Gaussian rule as we believe it represents a wider range of 
agent behavior. We propose to study such dynamics on 
both Erdos-Renyi (ER) [HI as well as scale-free (SF) [H 
networks. Thus a randomly selected node is seeded with 
a contagion of quality a and the subsequent dynamics is 
monitored. 

In order to take into account the heterogeneity induced 
by the presence of node with different connectivity, we 
consider the time evolution of the density of informed 
(infected) Ik(t) and uninformed (uninfected) U k (t) nodes 
of connectivity k at time t Q • These variables are related 
by means of the normalization condition: 



U k (t) + I k (t) = l 



(1) 



The total density of the informed nodes in a network 
is expressed by an average over the various connectivity 
classes; i.e., I{t) — ll k I k {t)p{k), where p(k) is the prob- 
ability of a node having k neighbors. These densities 
satisfy the following set of differential equations [TJ, [27| : 



dJ fc (t) 
dt 



kU k (t)f(k,a)e k (t) 



dU k (t) 
dt 



= -kU k (t)f(k,a)Q k (t) 



(2) 



(3) 



where the creation term is proportional to the node de- 
gree k, the probability U k {t) that a node with degree k is 
not informed, the local probability of acceptance (trans- 
mission) f(k,a) and the density Q k (t) of the informed 
neighbors of a node of degree k which received the con- 
tagion at time t. The last term is thus the average prob- 
ability that any given neighbor of a node of degree k is 
informed and received that contagion at the exact time 
t. For uncorrelated networks, the probability that each 
link of an uninformed node is connected to an informed 
node of degree kl is proportional to the fraction of links 
emanating from such nodes |2(]. By considering that at 
least one of the links of each informed node is connected 
to another informed node from which the contagion has 
been transmitted, one can show that e k (t) = 6(t) HE!) 
and in particular we obtain: 



Q(t) 



^ k ,(k' -l)p{k')I k ,(t)Q k ,(t) 
(k) 



Q k (t) = fcO(t-l) 



(4) 



(5) 



where (k) = Y> k ik'p(k r ) is the proper normalization factor 
dictated by the total number of links and Q k (t) is the 
probability that a node of degree k receives the contagion 
a at time t. 

A reaction rate equation for 0(t) can be obtained from 
Eqs. (3) and (5). In the initial spreading steps, one 
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neglects terms of order 1% and assuming, for early times, 
that the probability a node receives the contagion in its 
previous time step is one, leads to the following set of 
equations: 



dl k (t) 



dt 



d6(t) _ Q(t) 



kf(k,a)e(t) 



(G) 



dt 



(k) 



((k 2 f(k,a))-(kf(k,a))) (7) 



These equations can be solved with the uniform initial 
condition /fc(0) = Jo -C 1 yielding for the total average 
density of informed nodes I(t) = T^kp(k)Ik(t): 



I(t)=I [l + r(kf(k,a)}( 



W-1 
(k) 



)(e 



t/r _ 



where 



(k) 



(k 2 f(k,a))-(kf(k,a)) 



1)] (8) 



(9) 



is the growth time scale of an initial spreading in the 
network. 

Eq. (9) is a very interesting result since it highlights the 
importance of local dynamics in spreading through the 
appearance of /(fci,a) assumed to be a general function 
which could be probabilistic or deterministic. Perhaps 
more importantly, for a given class of local dynamics (a 
given /(fc, a)) the spreading rate depends on the quality 
of the contagion, a, being spread. 

We also note that our results is a generalization of the 
standard (mean-field) SI model for contagion spreading, 
i.e., for a global transmission probability, /(fc, a) — A, we 
obtain the standard result for heterogeneous mean-field 
models I3J: 



T = 



(k) 



(10) 



On the other hand, for a homogeneous network, where 
nearly most nodes have a degree close to the mean de- 
gree (k), our results can be simplified by considering 
(fc 2 ) = (k) + (k) |30|. Thus, the growth time scale in a 
homogeneous network can be obtained as: 



1 



(k)f((k),a) 



(11) 



which is to be compared with the usual r = 1/(A (k)) for 
homogeneous mean- field models Q. 



NUMERICAL RESULTS 

We now present our simulation results for two types of 
random networks, SF with p(k) oc fc -7 with 7 = 3 and 



an ER network with a Poissonian distribution. The net- 
work realizations used for the numerical simulations were 
constructed using the method introduced in [29[ in order 
to assure that no degree-degree correlations are present 
in any of the networks generated. We choose a = 0.1, 
and start by choosing a random site and infecting it with 
contagion a. The reported results are typically averaged 
over 1000 seeding and a separate (ensemble) averaging 
of at least 10 realizations of a given random network. 
Fig.l shows the time evolution of average density of in- 
fected/informed nodes I(t) on both networks with the 
same average degree (k) = 6 and network size N = 20000. 
Inset shows a log-linear plot of the early time behavior of 
I(t) indicating an exponential growth with a dependent 
time scales. Fig. 2 shows the growth time-scale r obtained 
from exponential fitting of numerical results along with 
the theoretical prediction (Eq.9), versus a. As can be 
seen, the numerical results recover the analytical calcula- 
tion with great accuracy. We see that the fastest spread- 
ing (smallest r) coincides with largest spreading (largest 
/). However, as a deviates significantly from the average 
quality of the network, (x) , r increases significantly, and 
thus greatly slowing down the spreading process. We also 
note that low quality contagion is more efficiently spread 
on SF networks, which can be expected, as SF networks 
have low average quality (x) — (k) /k max because of the 
arbitrary large k max . 

Based on the results shown in Fig.l, one might expect 
that initially fast spreading leads to a stationary state 
where a certain a-dependent percentage of the nodes are 
informed. However, this is not the case. The spreading 
actually never stops in our model, and due to random na- 
ture of our local rules, even the nodes whose quality do 
not match well with the incoming contagion can eventu- 
ally accept and thus transmit the contagion due to large 
number of exposures. This indicates that eventually all 
nodes are infected if one waits long enough. However, 
simulations show that the time required to reach such 
stationary state (I(oo) — 1.0) is very long and increases 
rapidly with system size N. What one typically observes 
after the fast exponential spreading is a quasi-stationary 
state where d//dt « on short time scales, along with 
a very slow increase to the eventual stationary state. Of 
interest, is therefore, the profile of such quasi-stationary 
states as a function of a. Fig. 3 shows such results for 
various times. We observe that the most efficient spread- 
ing occurs when the quality of the contagion matches the 
average quality of the network, i.e. a — (x). Signifi- 
cant spreading also occurs for a width of 2a about such 
maximum. For a SF network (Fig. 3(a)), average quality 
is nearly zero and the cutoff (at a = (x) + 2a) is rel- 
atively sharp. For ER network (Fig. 3(b)), the average 
quality can differ, but spreading falls off gradually with 
the width 2a about such an average. As different a's 
have different time scales, r(a), we also plot the profiles 
of quasi-stationary states after 100r(a). This is shown in 




FIG. 1: Averaged density of informed nodes versus time in 
(a) a SF network and (b) an ER network. The networks have 
the same size of N = 20000 and average degree (k) = 6. The 
insets show the log-linear plot of the same data for the first 
few time steps. Note that the average quality (x) is much 
smaller on the SF network. 
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FIG. 2: Measured time-scale r vs a as obtained from ex- 
ponential fitting in a SF network (circles) and an ER net- 
work (squares). Dashed lines indicate analytical calculation, 
Eq.(9). Parameters are the same as Fig.l 

Fig. 4 where the same type of behavior is seen as in Fig. 3. 

We next ask how the eventual stationary state is 
reached. We find that the quasi-stationary state evolves 
with slow dynamics I q s « t b ^ with time-dependent ex- 
ponents b(t) — J2 C jL2 c j/(^°S^y where c/s are a depen- 
dent. Fig. 5 shows the simulation results for various a's, 



FIG. 3: The density of informed nodes as a function of a for 
three different times on a (a) SF network and (b) ER network. 
Parameters are the same as Fig. 1. 
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FIG. 4: The density of informed nodes as a function of a at 
t — 100r(a) on a SF network and ER network. Parameters 
are the same as Fig. 1. 

along with a fit of the above function including the first 
three terms, fncluding higher order terms will make an 
increasingly better fit. Such multi-fractal slow dynamics 
to the eventual state is interesting in its own rights. The 
transition from fast initial exponential growth to the slow 
multifractal quasi-stationarity can be seen in Fig.l. 

As noted earlier, spreading slows down greatly for val- 
ues of a which significantly deviate from the average 
quality of the network (x) , because of exponential growth 
of r, see Fig. 2 and Eq.(ll). This leads to a low value for 
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FIG. 5: The long time behavior of the evolution of informed 
nodes on a (a) SF network of 7 = 3 and (b) ER network, both 
with the same mean degree (k) = 15. The dashed lines show 
the function t b ^ with b(t) = ^^^Cj/Qogty where Cj is a 
dependent fitting parameters and the fitting shows only the 
first three terms with R 2 « 0.98 while including higher order 
terms makes R 2 — y 1. 



FIG. 6: The size (N) dependence of the long time (t ~ 20r) 
value of the quasi-stationary state I q . s . = I(t « 20t) for 
(a) a SF network and (b) an ER network for contagion 
a = (x) + 0.25. The dashed line corresponds to I q . s . ~ 
with /3 = 0.57 for SF (a) and /3 = 0.13 for ER (b). The 
network parameters are as the same as in the Fig.l. 



Iq. s . which can be interpreted as low prevalence. How- 
ever, one can see that even if a does not significantly de- 
viate from (x), a finite, arbitrary r would lead, in a first 
few initial steps, to a finite number of infected/informed 
agents. It is therefore seen that low prevalence (i.e. 
small Iq.s.) could occur for large networks (N — > 00) for 
generic values of a. To show this we plot the value of 
I q .s.(t = 20r) for a = (ir) +0.25 (a = 0.1) for various sys- 
tem sizes. The results are shown in Fig. 6. Interestingly, 
for both network types one observes a scaling I q . s . ~ 
indicating low prevalence for large system sizes. We note 
that the exponent (3 is larger for SF networks than for ER 
networks showing a better tendency for low prevalence in 
SF networks. One might also wonder about the velocity 
dlq. s .(t)/dt as a function of N. Such results are shown 
in Fig. 7, where the velocity also decreases (to zero) on 
SF networks while it remains at a constant (albeit small) 
value for ER networks, again showing a better tendency 
for low prevalence on SF networks. 

As noted above, nodes remain active in this model and 
can eventually activate previously inactive neighbors due 
to repeated transmissions. This is due to the probabilis- 
tic nature of local rules considered here. One can there- 
fore wonder about the frequency of node's activity sim- 



ilar to how frequently users "retweef'or "comment" on a 
news in a social network [19,22]. We therefore define a 
node's activity as the number of times it has transmitted 
a given contagion, and plot the probability of such activ- 
ity for large networks over long times. Fig. 8 shows such 
results. Remarkably, for both types of networks and a 
wide range of contagion quality a (as long as the quality 
is not close to network average), a robust power-law be- 
havior is observed, P(a) ~ a~ s with 6 = 1.0. Note that if 
a approaches (x) then a deviation from pure power-law is 
observed with the tail becoming more heavily weighted, 
where eventually for a — (x) one observes a flat distribu- 
tion with a strong and sharp peak in the tail. Power-law 
behavior in agents activities has been reported in vari- 
ous empirical studies 20 23], with various exponents and 
various definitions of "activity". What seems to be im- 
portant in our model is the generic conditions for such 
power-law behavior, i.e. a wide range of a and network 
topology, as well as the seemingly exact exponent S = 1.0 
over many decades. 
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FIG. 7: The velocity of the quasi-stationary evolution 
(dj(t)/dt) at long time (t ~ 20r) versus system size N for 
(a) a SF network and (b) an ER network. The plot parame- 
ters are the same as in the Fig.l. 



CONCLUDING REMARKS 




FIG. 8: Distribution of node's activity for various contagion 
quality on both types of networks. The plots are shifted for 
clarity. Here, N = 20000 and (k) = 6 with 50 seeds and 50 
network realization for averaging. The dashed line is a line of 
slope —1 to help guide the eye. 



ing networks are heterogeneous in nature where (x) — > 
in such limits. Therefore, this result may explain the cu- 
rious popularity of tabloids, or the general prevalence of 
superstitious beliefs. 

Support of Shiraz University Research Council is 
kindly acknowledged. 



In this work we have proposed a general model of 
spreading phenomena which takes into account the qual- 
ity or strength of a given contagion. We also define a local 
probabilistic rule which models the fact that agents in- 
teract with contagions and transmit them based on their 
observed fitness. Our model reproduces some key fea- 
tures of spreading dynamics observed in empirical studies 
but absent in standard theoretical models. We find that 
the model exhibits quality-dependent fast initial spread- 
ing followed by a quasi-stationary state. The model also 
exhibits a robust power-law behavior of node's activi- 
ties. Both these effects have been recently observed in 
empirical studies of social networks 2(J 23 1, leading us 
to believe that such a model, despite its simplicity, is 
an accurate description of spreading dynamics on com- 
plex networks. Another important feature of our model 
is the observed low prevalence as a result of large system 
sizes. For arbitrary values of contagion quality a, and 
small <7, a finite number of agents are fit to receive the 
contagion and such fraction becomes ever smaller with 
increasing system sizes, more so on heterogeneous net- 
works. For large enough system sizes, only contagions 
which fit the average network quality a « (x) have any 
significant prevalence. We note that most naturally aris- 
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